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ABSTRACT 

We consider here the collective accretion of gas by globular clusters and dwarf 
galaxies moving through the interstellar medium. In the limit of high velocity 
and/or sound speed of the ISM, the collective potential of the cluster is insuffi- 

o 



cient to accrete significant amounts of gas, and stars within the systems accrete 
gas individually. We show, however, that when the sound speed or the relative 
velocity of the ambient medium is less than the central velocity dispersion of the 
cluster, it is accreted into the collective potential of the cluster prior to being ac- 
creted onto the individual stars within the cluster. The collective rate is strongly 
enhanced relative to the individual rates. This effect may potentially modify 
the white dwarf cooling sequence in globular clusters with low-inclination and 
^ ■ low-eccentricity Galactic orbits, and lead to the rejuvenation of some marginally 

surviving cores of globular clusters and nucleated dwarf galaxies near the Galac- 
tic center. Such effects will only occur rarely, but may explain the existence of 
clusters of young, massive stars near the Galactic center. 

Subject headings: galaxies: dwarf — galaxies: evolution — clusters: globular - 
clusters: open — hydrodynamics — methods: numerical 



1. Introduction 

Stars within individual open and globular clusters are coeval and have remarkably ho 
mogeneous chemical composition (see, eg. Penny & Dickens 1986; Richer & Fahlman 1986 
Stetson & Harris 1988; Ferraro et al. 1991, 1992; Kraft et al. 1992; Langer et al. 1992 
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Suntzeff 1993; Quillen 2002; Wilden et al. 2002). These properties provide strong con- 
straints upon theories of star formation and evolution within these objects. In the most 
metal deficient globular clusters such as M92, the total heavy element content is comparable 
to the output of a small number of supernova events. The observed upper limit on the color 
spread among their giant branch stars implies their metallicity dispersion AZ < 10 — 20%. 
Such small spreads place limits on the accretion of enriched interstellar medium by the main 
sequence progenitors, which can be strengthened further using differential spectroscopy. 

The accretion of the interstellar medium has more dramatic effects on both white dwarf 
and neutron stars. The white dwarf cooling sequences in both open and globular clusters 
are used as probes to determine the ages of these oldest stellar subsystems in the Milky 
Way galaxy (Richer et al. 1997, 1998; von Hippel & Gilmore 2000; Kalirai et al. 2001; 
Andreuzzi et al. 2002; Hansen et al. 2002). The intrinsic luminosity of the white dwarf is, 
however, typically < 1O~ 3 L . Accretion of ambient gas at a rate > 1O~ 14 M yr _1 would 
generate sufficient energy to modify both the color and magnitude of these degenerate dwarfs. 
Although the expected accretion rate of individual halo stars is several magnitudes below 
this value, it is not clear whether gas can be accreted first, at a sufficiently high rate, into 
the collective gravitational potential of individual clusters. The increased ambient density 
within the cluster would then enhance the accretion of the individual stars. 

The accretion of ambient gas by moving bodies is a classical problem. Many studies 
have been focused on the flow around compact stars with a point mass potential (cf Frank, 
King & Raine 2002). Although clusters have much larger masses than individual stars, their 
potential is relatively shallow. In this paper we consider the efficiency of accretion in these 
cluster potentials. We are interested in determining the accretion rate for clusters with a 
range of velocities relative to the ambient medium. Cases with high relative velocity and 
sound speed represent globular clusters belonging to the halo and thick disk populations. 
Cases with low relative velocity and sound speed are more relevant to young open clusters, 
globular clusters that have sedimented into the galactic disk, and the residual dense cores of 
globular clusters or nucleated dwarf galaxies near the Galactic center. 

In §2, we briefly recapitulate the relevant equations which describe the accretion flow. 
We deduce analytically a condition for the clusters to accrete gas collectively. Because the 
flow is multi-dimensional in nature, we adopt a numerical approach to verify our analytic 
approximation. We describe the numerical scheme and the range of model parameters and 
the results of these calculations in §3. These results are applied, in §4, to study the accretion 
of gas onto sedimented globular clusters as well as residual dense cores of clusters and 
nucleated dwarf galaxies near the Galactic center. In §5, we summarize our results and 
discuss their implications. 
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2. Accretion Onto a Shallow Potential 

The basic governing equations for the ambient gas are the continuity and momentum 
equations. In order to highlight the dominant physical effects, we neglect the energy equation 
and adopt a polytropic equation of state. In general, the accretion flow onto a gravitating 
body is known to be unstable and unsteady. Because we are interested in the time-averaged 
accretion rates for illustration purpose, we seek steady state solutions of the governing equa- 
tions such that 

V-pu = (1) 

pu ■ Vu = -VP - pV0 (2) 

where p, P, and u are the density, pressure, and velocity of the background gas. The 
gravitational potential of the cluster can be approximated with a Plummer potential, 

4>{R) = ~ GM \ n , (3) 
[Rl + R 2 ) 1/2 

where M d is the total gravitating mass, R c is the core radius, inside of which the density 
of the gravitating matter is approximately constant, and R is the distance from the center 
of the potential. The Plummer potential is relat ively centrally concentrated, and has been 



used to model globular clusters (IPlummer Ill91ll ). 



Equations (0Q) and (j2J) are mult i- dimensional in nature, and they are treated numerically 
in the next section. We can, however, obtain some approximate analytic solutions to delineate 
various regimes of interest. With a polytropic equation of state, in which P = Kp 1 where 
K and 7 are the adiabatic constant and the polytropic power index, we find from equations 
(ED and (ED that 

u • Vu = -c^Vlnp - £(1 + £ 2 )- 3 / 2 V£ (4) 

where 

i = R/Rc (5) 

is a dimensionless radius. The quantities u and c s are dimensionless speeds, defined as 

U= (GM d /R c y/*> (6) 

and 

° s = (GM d /R c y^ (7) 
where C s is the sound speed of the fluid. From the continuity equation (JTJ, 

u • lnp = Vu. (8) 
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We now consider three limiting cases. We first consider the possibility of very small R c such 
that both «(£ = oo) = <C 1 and c s (£ = oo) = <C 1. In this limit, the potential 
essentially reduces to that of a point mass and the accretion rate is given by Bondi solution 
(Bondi 1952; Hoyle & Lyttleton 1941; Shu 1992), such that 

M c = A( 7 )vrp(oo)(^ + M L) 1/2 J R| (9) 

where A is a constant of order unity, arid Uoo are, respectively, the non-normalized sound 
speed and flow velocity. The Bondi radius is given by 

GM d . 

R B = lr r> , , 2V ( 10 ) 



so that 

where = poo/^h and is the mass of the hydrogen atom. 

We note that the above form of Rb is derived for a point-mass potential, and so is only 
a good approximation for our models when Rb significantly exceeds R c . We now consider 
two limiting cases with finite values of R c . In case 2, we consider c s (oo) >> u(oo), such that 
the flow is nearly spherically symmetric and equations^ and (JHJ) reduce to 

1 U ^-7^V (12) 



din£ v« 2 -^/ V s (l + £ 2 ) 3/2 

The maximum value of £ 2 /(l + £ 2 ) 3 / 2 occurs at £ = £ m = \/2- Because c s is a monotonically 
decreasing function of £, the numerator of the right side of equation (j!2p remains positive 
definite provided 

C,(oo) > ^'\GM d /R c f' 2 . (13) 

Consequently, |tt| < c s throughout the flow. The ambient gas therefore collects into a quasi- 
static atmosphere, with a density distribution given by 



p(R) 



7 -l\ GM d 



if 7 y (R 2 + R 2 c y/ 2 



1/(7-1) 



1+, ~ 



^H(i+a i/2 



1/(7^1) 



(14) 



where is the density of the ambient medium. If = C s (oo) > (GMd/Rc) 1 ^ 2 , the density 
enhancement throughout the cluster is very small. 

The relative motion of typical globular clusters with respect to the interstellar medium 
is generally supersonic, i.e. <C u^. The flow pattern for passage around the cluster with 
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a large relative velocity, u rXl} is multi-dimensional and complex. In line with the conventional 
treatment of Bondi-Hoyle accretion flow, it is intuitively tempting to modify the transition 
condition to be 

C 2 00 + u 2 00 >GM d /R c . (15) 

When this condition is satisfied, it is more appropriate to use the Bondi-Hoyle accretion 
formula for individual stars 

M. = B(~,)irp(K)(GM,nCi + ul)-W 

= 10»B(7) (j^fc*) (fc)' (l + f)' 3 ' 2 (ufT^r S «"» <16) 

where -6(7) is a constant of order unity (it is not the same as A due to the difference in the 
potential) and M* is the mass of the star. The total accretion rate of the cluster is 

M t = NM* (17) 

where the number of stars in the cluster is N ~ M^/M*. But when equation ([15]) is not 
satisfied, M c in equation is the more appropriate rate for the cluster. Note that 

M c ~ NM t » M t . (18) 

This inequality clearly demonstrate that the collective effect of the cluster would be greater 
than sum of the individual stellar contributions if the it moves subsonically through the 
interstellar medium. 

In the next section, we will carry out 3-D numerical computations to verify the condition 
in equation (Tl5|) . Here we provide some analytic approximation for the limit that » C^. 
Without the loss of generality, we consider a case 3 in which the ambient gas is approaching a 
cluster in the x direction and the downstream flow velocities can be written as u x = + 5u x 
and u y = 5u y . Neglecting contribution from the pressure term in equation P|, we find that 
along the streamline with an initial impact distant Ri, 

8ux GM d 

In the absence of shock dissipation, the potential vorticity, u/p, and the Bernoulli energy, 
\u\ 2 /2 + W + 0, are constant of motion along stream lines. In the limit that the enthalpy 
W = J dPj p is small compared with the kinetic and potential energy, the Bernoulli constant 
implies that ^Uy/u^ ~ \5u x /uca\. The deflection angle of the stream line Sy u /uoo is small 
in the limit u 2 ^ > GMd/R c - We also note, from the conservation of the potential vorticity, 
that, in the high velocity limit, the modification of p is small despite the slight convergence 
of the streamlines. 
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3. Numerical Models 
3.1. Method 

The above analytic approximation is greatly simplified. The flow is actually multidi- 
mensional in nature. We therefore proceed by considering multidimensional simulations of 
accretion by a defined potential moving relative to a dense gas cloud. The results found here 
shall be used to validate the analytic estimates of § 2. 

The numerical code we use is Cosmos. It is a massively parallel, multidimensional, 
radiation-chemo-hydrodynamics finite difference code developed at Lawr ence Livermore Na- 



tional Labo ratory. This scheme, and tes t s of t he code are described in lAnninos fc Fragile 



( 120031 ) , and lAnninos. Fragile fc Murray I (120031 ) . It has also been used to st udy the evolu- 



e et al. 
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Murray. Dong fc Lin I (120031 ) . 



tion of supernova-enriched material in dwarf galaxies by iFragi 
Roche Lobe overflow from cluster and galactic potentials by 
The reader is referred to those papers for more details on Cosmos. We discuss here the 
settings used for the current work. 

For a spherically symmetric potential, the unperturbed flow would be expected to be 
two-dimensional. In flows where Woo ~ Coo, however, an axisymmetric flow pattern may 
be unstable (Landau & Lifshitz 1959; Batchelor 1967). Our objective is to determine the 
time averaged accretion rate rather than to identify any flow instability. Because Cosmos 
is written in Cartesian coordinates, however, we cannot impose cylindrical symmetry upon 
the problem, and so the simulations are carried out in three-dimensions. So as to reduce 
computational expense, the models are quadrants, with one-fourth of the regions included. 

We wish to examine the ability of a potential to accrete gas as a function of the relative 
speed of the potential through the gas, and the gas temperature. To improve the controlled 
nature of the models, we do not include radiative heating or cooling. The gas, instead, 
evolves adiabatically. The effects of radiative equilibrium are approximated by having the 
gas evolve with an adiabatic constant 7 = 1.01, giving approximately isothermal behavior. 
In cases where sufficient gas is accreted for it to become self-shielded, cooling could decrease 
the temperature of the gas significantly, potentially enhancing the accretion rate beyond the 
values computed here. 

Rather than examine the motion of a potential through a static cloud, we take the 
simpler, equivalent approach of viewing the encounter within the frame of the cluster, in 
which a dense cloud of gas sweeps over the fixed potential. The model clouds are extremely 
simple. They are slabs, with thickness I = 100 pc, moving across the problem with speeds 

Of Woo- 
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As discussed in § 2, the important quantities are the relative speed of the cloud and the 
potential, u^, and the sound speed of the gas, C s , and these quantities are varied between 
the models. The gas within the dense slab has temperature selected to give the desired value 
of C s and a density chosen to cause it to be in pressure equilibrium with a hot background 
gas. The hot background gas in the problem has T = 10 6 K and p = 1.6 x 10~ 25 g cm -3 . 
The pressure of the gas, nT = 10 5 , is chosen to be comparable to that expected within the 
central regions of galaxies, or in dense cooling flows. 

As stated above, accretion by a potential moving relative to a gas cloud is complex. In 
the frame of the potential, the gas streamlines are bent towards the cluster center. Some 
shall intersect the center, while others converge along a line behind it. The convergence 
speed of the gas determines the reduction in its velocity relative to the potential due to 
shocks, and therefore whether or not the gas is accreted. In order to ensure that we are 
accurately measuring the mass of accreted gas, we carry out our models until the dense slab 
of gas has completely swept past the potential. The increase in the amount of gas within 
the potential radius as compared to that present initially, M acc , is computed, and listed for 
each of our models. If the value of M acc is seen to be changing at the end of the simulation, 
the model is run for additional time. The minimum time for the simulations is a few times 
i/«oo«100Myr (^). 

For the models examined here, the total gravitating mass, M d = 7 x 10 5 M , and the 
core radius of the potential, R c = 10 pc. To approximate the effects of external gravitational 
fields, the potential is flattened (the gravitational force goes to zero) beyond an assumed 
tidal radius Rt = 80 pc. The central velocity dispersion a = (GM^/ 'i? c ) 1//2 = 16 km s _1 , 
similar to the values within many globular clusters (Pryor & Meylan 1993; Djorgovski 1993). 
In some centrally condensed clusters such as M15, a can reach 20 km s _1 . In some nucleated 
dwarf galaxies, a > 30 km s~ x (Peterson & Caldwell 1993; Geha, Guhathakurta, & van 
der Marel 2002). In contrast, a in the Pleiades and Hyades clusters are 0.6 and 0.2 — 
0.4 km s" 1 , respectively (Perryman et al. 1998; Chen & Zhao 1999; Madsen, Dravins, & 
Lindegren 2002). Despite these differences, the discussion in the previous section indicate 
that the governing equations can be normalized with dimensionless parameters and the 
results obtained therefrom can be scaled to the appropriate limits. 

Self-gravity of the gas is not included in the models. This approximation should be 
adequate for most of our models, for which the accreted mass is less than the mass responsible 
for the potential. For two of the models discussed below, the accreted mass exceeds that 
of the potential. For those models, the inclusion of self-gravity would further enhance the 
accretion rate. Both because of this, and the additional cooling discussed above, the accreted 
masses for those two models should be treated as lower limits. 
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The physical dimensions of the models are 400x200x200 pc, with resolutions of 2 pc. 
The core radius, R c , is therefore resolved by only 5 zones. While not highly resolved, the 
fraction of the accreted mass contained within the core is generally small, and so the limited 
resolution within the cores does not affect our conclusions. 



3.2. Model Results 

The parameters of the models which we have run are shown in Table [TJ In the table 
are listed, for each model, the speed of the gas relative to the potential, Uoo, the sound 
speed of the gas, C s , the Bondi radius (as defined in equation [TD]), R B , the crossing time of 
the potential across the dense slab at Mqo, t c = Z/ttoo, the sound crossing time across Rb, 
t s = Rb/C s , the amount of gas predicted by approximate theory to be accreted (see below), 
M p , and the actual amount of gas accreted in the model, M acc . The value of M acc is taken 
as the increase in mass contained within Rb relative to that present at the beginning of the 
simulation, and we ensure that the value is not changing by the end of the simulation (see 
above). 

The parameters of the models are chosen so as to span the range of possible behaviors. 
In Model 1, > a, while in Model 6, C s > a. In both models, little accretion is expected, 
but in Model 1 this is because of the high relative velocity, while in Model 6 it is due to the 
high sound speed of the gas. In Model 2, ~ o » C s , while in Model 5, C s ~ o » Uoo. 
Both models represent marginal conditions for accretion, but again for different reasons. In 
Model 3, Moo = C s < a, and so that model is expected to represent an ideal situation for 
substantial accretion. In Model 4, = C s ~ a, and so, again, that model represents a 
marginal situation for accretion. 

Snapshots showing the evolution of Model 2 are shown in Figure HJ In the figure, the 
dense slab of gas can be seen sweeping over the gravitational potential, which is visible due 
to its effect upon the density of the gas moving past it. Because ~ 05, some gas is pulled 
ahead of the cloud by the potential, and a large fraction of the gas within the tidal radius 
of the potential is accreted during the cloud passage. Once the dense slab of gas is past 
the potential, a small, dense core of accreted gas remains in the center of the potential. In 
the models which accrete significant amounts of gas, therefore, the gas rapidly evolves to a 
condensed configuration, such as is likely to lead to cooling of the gas, and greatly enhanced 
accretion by the individual stars within the system. 

We may make a crude analytical estimate of the accreted mass if we assume that the 
potential acts as a "cookie cutter," and that gas closer to the center of the potential than 
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Rb is accreted during the passage of the potential through the dense slab. The amount of 
gas that would be predicted to be accreted is then given by 

M p = nlR 2 B p s , (20) 

where p s is the density of the gas in the slab. The calculated values of M p are shown in 
Table [U 

As can be seen from the table, the values of M acc and M p are in fair agreement for 
some models, but not others. As indicated in §2, the expression of Rb in equation (TTDT) 
is derived for a point-mass potential, and so is only a good approximation for our models 
when Rb significantly exceeds R c , a condition violated in Models 1 and 6. For those models, 
the Bondi radius computed assuming a point mass potential lies within the core, where 
the Plummer potential actually flattens out. The true potential therefore never becomes as 
deep as required for substantial accretion, and M acc <C M p . Because Uoo > C s in Model 1, 
equation ffTTl) is more appropriate to represent the gas accretion. We could not, however, 
include the fine grain structure of the potential due to individual stars in our simulations. 

In the case of Model 3, Rb > Rt- Beyond R t , the gravity of the potential can have no 
direct effect upon the gas, and so Equation [20] represents an overestimate of the amount of 
gas that is expected to be accreted. Had we used R t instead of Rb in calculating the fiducial 
value of M p , we would have found M p = 5.6 x 10 7 M Q , much closer to the value found for 
M acc . 

For the remaining models, M acc > M v . That discrepancy is likely to be the result of the 
simplistic assumption made in deriving M p , that no hydrodynamic motion beyond R t results 
from the accretion. That assumption is equivalent to assuming that t s > r c . In fact, for most 
of the models, t s is at least comparable to, and often substantially smaller than r c , indicating 
that the gas within the dense slab shall have time to undergo significant adjustments to a 
new hydrostatic equilibrium in response to the accretion. As gas is accreted by the potential, 
therefore, the resulting pressure gradient drives gas inward from beyond R t , increasing the 
effective radius from which gas may be accreted. The resulting enhancement in the accretion 
causes M acc > M p , even for Model 6, for which Rb < R c - 

The results of these models indicate that the results of Section [2] represent reasonable es- 
timates of the accretion rates onto cluster potentials, and that the ratio (Cf + u 2 ^) / {GM d / R) 
is crucial to determining the importance of accretion by the global potential. The greatest 
differences between the analytic and numerical results occur at low values of M, where the 
accretion is insignificant. 
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4. Astrophysical Applications 

We now consider three scenarios in hich accretion may affect the evolution of stellar 
clusters. These include potential affects upon the metallicity dispersion and white dwarf 
cooling sequences of globular clusters, the stellar populations of clusters in the Galactic disk, 
and the rejuvenation of old clusters near the Galactic center. 

4.1. Metallicity Dispersion and White Dwarf Sequence in Globular Clusters 

Typical globular clusters in the Galactic halo have relatively high orbital inclination 
and eccentricity. In contrast, the interstellar medium consists of either warm (10 2 K) atomic 
or cold (10 K) molecular gas on nearly coplanar, circular Galactic orbits. The relative 
velocity of the halo clusters to the interstellar medium is in the range of 100 km s _1 . For 
these clusters, the collective effect of their potential is clearly insignificant, and cluster stars 
accrete as individual entities. 

Scaling with an average number density of 10 2 cm -3 for atomic clouds, we find from 
equation HTB]) that M* ~ 10 9 g s _1 . Over the age of the Galaxy, the total amount of mass 
accreted is ~ 10~ 6 M Q . Even for the most metal-deficient ([Fe/H] ~ —4) stars known, the 
acquisition of this amount of interstellar gas, with a solar metallicity, would modify the 
metallicity of the star by less than 1%. Variations in Lithium abundance has been observed 
among subgiants in halo globular clusters (Deliyannis, Boesgaard, & King 1995; Castilho 
et al. 2000). Main sequence evolution almost certainly has depleted all the Lithium in the 
envelopes of globular clusters' main sequence stars and the above estimate for the accretion 
rate is too small to account for the observed data. But, Lithium is produced in the core 
through hot bottom burning and dredged up to the surface of stars on the asymptotic and 
red giant branches (Cameron & Fowler 1971; Sackmann, Smith, & Despain 1974; Sackmann 
& Boothroyd 1992). The consumption of a Jupiter-mass gas planet has also been suggested 
as a mechanism to replenish the surface Lithium content (Laws & Gonzalez 2001; Sandquist 
et al. 2002), though the depletion time scale for Lithium is relatively short. 

For white dwarf stars, accretion at the above rate would lead to the release of energy at 
a rate L a ~ 1O~ 7 L , which is much smaller than that expected from the cooling of the oldest 
Galactic white dwarfs (Hansen 1999). Thus the accretion from the interstellar medium needs 
not be taken into account in the determination of the white dwarf cooling sequence. 

In order, therefore, for accretion to have a significant effect upon either the metallicity 
of cluster stars, or upon the the color and luminosity of white dwarfs, Mqo must be reduced to 
~ 10 km s^ 1 , in which limit the effects of collective accretion by the overall cluster potential 
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would become important. 



4.2. Sedimentation of Clusters into Galactic Disks 

Although for typical halo clusters is too large for them to accrete any significant 
amount of matter, their orbital properties may evolve. As halo clusters pass through the 
Galactic disk with a speed u g , the wake of disk stars induce a drag force (dynamical friction), 
given by 

F d = AC d 7cG 2 Mjp d u- 2 (21) 
where p d is the stellar density in the galactic disk, the drag coefficient, 



C, = InA 



erf(X)-^e-* 2 

7T i /2 



(22) 



where A = b m3iX /R c , X = u g / (\/2er), and 6 max is the maximum effective impact parameter, 
which can be taken to be the density scale height H d of the Galactic disk (Chandrasekhar 
1943; Binney & Tremaine 1987). Because gas and stars are concentrated near the midplane, 
the rate of dissipation of the cluster's orbital energy, averaged over its Galactic orbital period, 



Pg, IS 

'-'(SO 

where E = u 2 g /2 is the kinetic energy relative to the disk and u z is the velocity normal to 
the galactic disk. Clusters which cross the Galactic disk with large velocities would sediment 
into low inclination and low eccentricity orbits on a time scale u 2 /4E (Artymowicz, Lin & 
Wampler 1993). Over a time interval, AT, the critical condition for clusters to attain small 
u g and u z , and therefore Uoo is 

^V^<32C>G 2 ^%^ (24) 
VJ V c M g {af Pg 1 ' 

where the Galactic circular velocity V c = [GM g (a) / a] at a radius a is determined by 
the Galactic mass M g (a) contained within it, and = p d H d is the surface density of the 
disk stars. For V c = 220 km s _1 at a = 10 kpc, equation (|24"I) implies that clusters with 
Ug ~ u z < 0.3V C or Ug = V c and u z < 10~ 2 V C would sediment into the disk within the life 
span of the Galaxy. Because AT is a decreasing function of u z and u g , the inclination and 
eccentricity of the clusters' Galactic orbits quickly vanish, leading to a small u^. 

Based on the present phase-space distribution of globular clusters, several authors (Fall 
& Rees 1977; Keenan 1979; Ostriker & Gnedin 1997; Vesperini, 1998; van den Bosch et al. 
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1999; Hideaki & Yoshiaki 2003) have suggested that dynamical friction may have caused 
clusters with small initial orbital radii to undergo significant orbital evolution. There also 
exist clusters in a thick disk population, which have smaller values of than typical halo 
clusters. 

For clusters that are not on orbits highly inclined to the disk, collective accretion may 
become a significant factor in their evolution. As is reduced to ~ 10 km s _1 , M* 
approaches to 10 13 g s _1 . More importantly, < (GM d /i? c ) 1//2 , and M c approaches 10 25 g 
s _1 . The gas density p c within the clusters' potential would increase until 

M* ~ (M d /M*)M C ~ 10- 7 {M*/M Q ) 2 M Q yr" 1 . (25) 

Because clusters have nearly isothermal stellar distribution function, Uoo relative to the gas 
in the cluster is comparable to a, and p c ~ (M^/ M*)p g ~ 10~ 16 g cm -1 . Gas within these 
clusters would achieve high densities. 

Within this environment, the cluster stars would accrete gas and gain significant mass. 
Within one Galactic orbit, the observed metallicity of the stars could attain the solar value 
present in the interstellar medium. Gas accumulating within the core of the cluster potential 
may also trigger new star formation. With high temperatures and densities, the initial mass 
function may be heavily skewed toward the high- mass range (McKee & Tan 2003). As 
star formation proceeds, negative feedback from processes such as photoionization, stellar 
winds, and sup ernovae can modify the gas concentration, and terminate star formation 



(IDong. Lin, fc Murray 1 120031 ) . These negative feedback effects cannot, however, prevent re- 
accumulation of gas around the sedimented clusters after the newly formed massive stars 
evolve off the main sequence. The distinguishing properties of these sedimented clusters 
from open clusters are their centralized cores and rich population of white dwarfs, despite 
their gas rich environment and the pre-main-sequence color magnitude diagram. 

There are no stellar clusters observed in the Galactic disk which bear these anticipated 
properties, as would be expected given the relatively small region of phase space occupied 
by clusters whose orbits would be heavily modified by dynamical friction. At early epochs, 
however, many more such clusters would be expected, as a natural extrapolation of the cold 
dark matter scenario in which merger events and the subsequent dynamical friction would 
lead to the dynamical evolution of not only globular clusters but also the cores of dwarf 
galaxies (Moore et al. 1999; Klypin et al. 1999; Metcalf 2002; Metcalf & Zhao 2002; Murray 
et al. 2003). Star formation via accretion onto stellar and dark matter potentials may 
therefore have played an important role in the early evolution of galaxies. 

Accretion at lower levels than the maximum estimated above would still affect the white 
dwarf cooling sequences of clusters. Clusters with orbits on small radii, and close to the 
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galactic plane might therefore be expected to have white dwarf cooling sequences brighter 
than expected for their ages. 

4.3. Accretion Onto Remnant Cores of Globular Clusters and Nucleated 

Dwarf Galaxies 

Near the Galactic center, several young massive star clusters have been recently discov- 
ered (Nagata et al. 1990, 1995; Okuda et al. 1990; Figer et al. 1999). Besides their age, 
these cluster are unusual in that 1) the initial mass function of stars within them is highly 
skewed to the high mass range (Figer et al. 2002; Stolte et al. 2002); and 2) they have 
short dynamical life expectancy (Kim, Morris, & Lee 1999; Kim et al. 2000; Portegies Zwart 
McMillan, & Gerhard 2003; MacMillan k Portegies Zwart 2003). One cluster, the Arches 
Cluster, contains up to 7xl0 4 M within a radius of 0.23 pc, and has a velocity dispersion 
of up to 22 km s^ 1 (Figer et al. 2002). In addition, the Quintuplet Cluster, the cluster 
around the Galactic center itself, and several other candidate clusters (Law & Yusef-Zadeh 
2004) contain x-ray point sources. The presence of clusters of young, short-lived stars near 
the Galactic center is remarkable, given that it is an environment that is extremely hostile 
to the formation of star clusters. 

Assuming that the above objects are all young and short-lived clusters, and that there is 
nothing unusual about the current epoch, then more than 10 4 such clusters would be expected 
to have existed over in the lifetime of the Galaxy. The inferred heavy element generation by 
the prolific production of massive stars would enrich the Galactic center region to well above 
the present level, unless gas is effectively removed from that region (cf Figer et al. 2004). 

Based on the results of the current work, we suggest an alternative scenario for the 
origin of these extraordinary clusters. Under the influence of dynamical friction, globular 
clusters undergo orbital decay towards the Galactic center. As discussed above, dynamical 
friction is not important to the evolution of most clusters. It is computed, however, to have a 
significant effect for clusters within a few kpc of the Galactic center (e.g. Fall & Rees 1977). 
The evolution of such clusters, and accretion of gas by them, may explain the existence of 
young clusters of massive stars near the Galactic center. 

During the course of the clusters' orbital decay, the Galactic tidal field gradually in- 
creases and their tidal radius gradually shrinks. Stars close to or outside of the tidal radius 
become detached while those near the dense core are retained (Oh & Lin 1992; Oh, Lin, & 
Aarseth 1995). Stellar loss also shortens the two-body relaxation time scale, enhances the 
central concentration, and shortens the tidal disruption time scale (Kim & Morris 2003). 
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With a core density comparable to that of post-collapse clusters M15 (Guhathakurta et al. 
1996) and M30 (Yanny et al. 1994), systems with more than 10 4 stars can persist and be 
tidally preserved near the present location of Galactic center clusters. Less centrally dense 
clusters such as M13 (Cohen et al. 1997) are likely to be disrupted by the Galactic tide 
outside 30 pc. 

We suggest that Arches and Quintuplet clusters are the rejuvenated cores of old globular 
clusters, which have retained their integrity as their orbits decay to the vicinity of the Galactic 
center. Of those, however, only a fraction have velocities relative to the atomic and molecular 
clouds that are sufficiently small for collective accretion to become significant. The velocity 
dispersion of the stars and gas clouds within 30 pc from the Galactic center is oq ~ 50km 
s _1 (Genzel, Hollenbach & Townes 1994). The central velocity dispersion of pre-collapse 
cores is typically a* ~ 10km. The fraction of clusters with + u 2 ^ < al (cf the condition 
set in equation [T5l) is ~ erf^/ag) ~ v1/&g> which is on the order of a few percent. 

For clusters having Md ~ 10 4 M Q and sufficiently small speed relative to the ISM, we find 
from equation (TiTj) that the gas accretion rate onto the cluster M c > 10 23 gm s _1 in regions 
where the density of the ambient gas near the Galactic center is > 10 4 cm -3 (Launhardt, 
Zylka, & Mezger 2002). The accumulation of gas near the center of the cluster potential 
enables solar-type stars to accrete at a rate ~ 1O~ 6 M yr _1 . Because the accretion rate is 
proportional to M 2 , the growth rate of the stars accelerates rapidly as they accrete mass. 
That rapid acceleratoin allows stars with masses of 20 M to be built up within ~1 Myr, by 
accretion onto stars with initial masses comparable to that of the Sun, a timescale shorter 
than the lifetimes of the massive stars formed via accretoin. With such a small velocity 
relative to the dense gas, the clusters do not drift more than a few pc from the orbits of any 
clouds within 1 Myr, which would allow the low-mass mature stars to be rejuvenated into 
young, massive stars. 

The main advantage of this scenario is that it does not require the rapid formation 
of clusters with unique initial mass functions. In comparison with star forming regions 
in the Solar neighborhood, the Galactic center is a challenging environment for efficient 
star formation. Our scenario naturally bypasses the stringent conditions for forming the 
progenitor clouds in regions where the external tidal effect is strong and heating is intense. 
This scenario requires a combination of clusters surviving tidal destruction as their orbits 
decay, and encountering dense gas with sufficiently low relative velocity to lead to significant 
accretion. Such a combination of events is certainly relatively unlikely. However, in the cold 
dark matter scenario, the initial population of clusters and dwarf galaxy cores near the center 
of the Galaxy was much larger than today, such that the population of faint, old clusters near 
the Galactic center would be significantly larger than predicted from the current population 
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of globular clusters. A fairly unlikely sequence of events is therefore needed if we are to avoid 
having many more systems such as the Arches and Quintuplet clusters today. 

Additionally, the centers of dense clusters are populated with both white dwarfs and 
neutron stars (Yanny et al. 1997). Accretion onto these dense stellar cores at the inferred rate 
can lead to the onset of luminous x-ray sources. Such a process may account for the distinct 
point sources of X-ray emission found by Chandra in these clusters (Laws & Yusef-Zadeh 
2004). 

5. Summary and Discussion 

In this paper, we examine in detail the accretion of ambient gas by a stellar system. 
Although the gravitational potential of the individual stars are point mass, that of the cluster 
is softer. We showed that if either the relative speed of the cluster or the sound speed of 
the ambient gas is large compared with the velocity dispersion of the core, the accretion 
by the cluster is inefficient. Individual stars accrete gas as though they move through the 
interstellar medium independently. But if the relative speed of the cluster and the sound 
speed of the ambient gas is less than that of the internal velocity dispersion, gas is accreted 
into the cluster potential collectively and rapidly. Accretion by individual stars is then 
enhanced greatly relative to their rate of accretion directly from the ambient gas. 

We suggest that this process may be important in inducing chemical inhomogeneity 
and modifying the white dwarf cooling sequence some globular clusters whose orbits happen 
to lie close to the Galactic plane. We also speculate that this process may have caused 
some tidally stripped cores of globular clusters near the Galactic center to accrete gas and 
to rejuvenate their member stars, resulting in the formation of clusters of young, hot stars, 
such as are seen near the center of the Galaxy. 

This work was performed under the auspices of the U.S. Department of Energy by 
University of California, Lawrence Livermore National Laboratory under Contract W-7405- 
Eng-48. This work is partially supported by NASA through an astrophysical theory grant 
NAG5-12151. 
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Table 1. Models 
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Fig. 1. — The density evolution of Model 2. The model is shown at the times 0, 4.5, 6.0, 
7.5, 22.5, and 30 Myr. The dense cloud, modelled as a simple slab of gas, sweeps from left 
to right across the accreting potential, centered in the model. Because of the relatively deep 
potential, and slow speed relative to the cloud, the potential has a significant effect upon the 
dense gas. Some of the gas is initially pulled ahead of the rest of the cloud into the potential, 
and much of the gas within the tidal radius is accreted as the cloud passes by the potential, 
where it forms the small, dense core seen in the final frame. 



